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Abstract 

In this paper, we study four-body systems consisting of diquark-antidiquark, 
and we analyze diquark-antidiquark in the framework of a two-body (pseudo-point) 
problem. We solve Lippman-Schwinger equation numerically for charm diquark- 
antidiquark systems and find the eigenvalues to calculate the binding energies and 
masses of heavy tetraquarks with hidden charms. Our results are in good agreement 
with theoretical and experimental data. 
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1 Introduction 


The idea of existence of tetraquark hadrons (two quarks and two antiquarks) was initially 
raised about twenty years ago by Jaffe. He utilized one of the initial versions of the bag 
model to study tetraquark spectroscopy of q 2 q 2 in which q was a quark lighter than charm 
quark. The MIT bag study revealed a dense spectrum of tetraquark states in the light 
sector pp. Later, tetraquark systems were examined in potential models and flux tube 
models [2j. The notion of diquark (two-quark system) is of use in describing the hadron 
structure and particle interactions at high energies. 

According to the quark model, hadrons are made up of quarks. Mesons consist of a 
quark and an antiquark in a bound state, such as light scalar mesons and some charmed 
mesons [3] and baryons are composed of three quarks in a bound state. Their structure 
was also shown to contain diquarks EaElE]. 

Based on the diquark model, two quarks join and build a color-anti-triplet in a bound 
state. As a simplified image, a diquark is viewed as a point particle having the quantum 
number of two quarks. More generally, a diquark is a system of two quarks considered as 
a pair. For example, a two-quark correlation in a hadron with at least two quarks will be 
a diquark system. In its ground state, a diquark has positive parity and may be an axial 
or a scalar vector. 

Gell-Mann [7] originally proposed the existence of diquarks. Based on the fundamental 
theory, the concept of diquark was developed, and it was even used to account for some 
experimental phenomena [8]. Ref. [9] explored the exotic state X(3872) via QCD Sum 
Rules. This study treated the hadronic state as a hidden-charm-tetraquark state with the 
quantum number I G (J PC ) = 0 + (l ++ ). Chen and Zhu [10j used the same interpolating 
current to investigate the l + ~ tetraquark state and found its mass to be (4.02 0.09)GeV. 

Maiani et al.|3j studied the masses of hidden-charm diquark-antidiquark systems tak¬ 
ing into account the constituent diquarks’ masses and spin-spin interactions, but Ebert et 
al. [I] employed the relativistic quark model based on the quasi-potential approach in order 
to find the mass spectra of hidden-charm tetraquark systems. Unlike Maiani et ah, they 
ignored the spin-spin interactions inside diquark and inside anti-diquark. We, however, 
considered tetraquarks as two-body systems, and spin-spin interactions were ignored. We 
used the diquark-antidiquark picture to reduce a complicated four-body problem to two 
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simpler two-body problems. The paper is organized as follows: 

In the first part, the bound states of four-quark systems are investigated in the framework 
of a pseudo-point two-body system. We explain Gauss-Legendre method in the second 
part. In the third part, we calculate the binding energy of heavy tetraquarks with hidden 
charms, and the last part is devoted to conclusions. 


2 Tetraquarks represented through two-body prob¬ 
lems 

Among the tetraquark states, those consisting of diquark-antidiquark are of interest to this 
study. To describe tetraquarks, we have taken diquark-antidiquark as if the pseudo-point 
diquark was a two-body system. Such interpretation helps reduce a complex relativistic 
problem to a simple two-body problem (Fig.l). 


Figure 1: Tetraquark systems 


Due to the wide variety of heavy diquarks, we have narrowed the study down to hidden- 
charm diquarks. These particles consist of a charm quark and a light quark ( u,d,s ). In 
order to use unrelativistic Schrodinger equation with tetraquark systems, we take heavy 
diquark and solve homogeneous Lippman-Schwinger equation numerically [Tlj |T2jJ (which 
is the integral form of Schrodinger equation) for two-body systems composed of diquark- 
antidiquark . 

Schrodinger equation for a two-body bound state with the potential V (which is as- 
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sumed to be energy-independent) runs as the following integral equation [II] : 


0b >— GqV | 0ft > 


( 1 ) 


Go is the propagator of a free particle. In configuration space, Eq.([I]) turns out as: 

] ro o 


0 6 (r) = -m— [°°dr'r' 2 [ dx’ [ * dcj) 1 
4ttJo J -i Jo 


r i r 2 -rr exp(—J m\E b \\r — r'l) 
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where E b stands for the binding energy of the two-body bound system (diquark+antidiquark). 
The wave function will be: 


0b(r) = j dr J dx M(r,r , x )0&(r ) 


where: 


1 exp((-Jm\E b \)Vr 2 + r 2 - 2rrf) 
M(r,r ,x) = —znm— - v - r - V(r ) 


4n \/r 2 + r 2 — 2 rrx 

Eq.(l4ll is of the following eigenvalue form: 


( 3 ) 


( 4 ) 


K(E b ) |0 6 >= \(E b )\ii> b > 


( 5 ) 


To determine the binding energy, we diagonalize the kernel in our calculations to obtain an 
eigenvalue A = 1, which is indicative of a physical system. The masses of the constituent 
quarks are as follows: 


rn s = 0.5 GeV m u = m ( i = 0.33GeE m c = 1.55GeE (6) 


The other required data include an arbitrary potential, r-cutoff, and reduced mass of the 
diquark-antidiquark. .R-cutoff is supposed to be the point at which the potential tends to 
zero. 


3 Gauss-Legendre method 

The eigenvalue equation (5) is solved through iteration method (direct method) p3jj- To 
discretize the integrals, Gauss-Legendre method ra is employed. Gauss lattice points 
for r,rf,xf are supposed to be 100. In Gauss-Legendre method, each integral of [-1,+1] 
interval is treated as: 

/■+ 1 

/ f(x)dx = J2 w if( x i) ( 7 ) 

J ~ ] i=i 
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where Xi denotes the roots of the type-one order-N Legendre function, and Wi are the 
functions of point weight. The following variable change is used to transfer the integration 
interval of r/ from [0 ,r max \ to [—1, +1]. If the integrals are discretized, then: 


1 + x 


r = r r 


NL NL 


Eq.( 


M r ) = - 2 ™f/2EEK^/. 

j= i*=i 

could be reformulated as: 


, , , 2 e xp(—yJm\E b \p(r, r', x')) 


p ( r , 




( 8 ) 

(9) 


N r' 

M r ) = N ( r > r i)^£»( r i) ( 10 ) 

2=1 

where: 

^ 2 exp(-Y /, m|E f) |p( 

p(r, 

Matrix N is diagonized to find A = 1 in the eigenvalue spectrum. The energy correspond¬ 
ing to A = 1 will be the system’s binding energy. 


iV(r.r') = 


r, r', x')) 


V(r[) 


( 11 ) 


4 Binding energy and masses of heavy tetraquarks 
with hidden charms 

In this part, we use a potential presented in | [T5] . In this potential, diquark interaction 
results from gluon field (spin-spin interaction in the potential is ignored). The potential 
is of this form: 


V (r 

) Vcoulir) + V con f(v^j 

(12) 

Vcoul 

4 F 1 (r)F 2 (r) 

(13) 

( r ) - o Ola 


3 r 


Vconf{r) = Ar + B 

(14) 


where F{r) is the form factor, which enters the vertex of the diquark-gluon interaction. 

F(r) = 1 - e« r - Cr2 (15) 

A = 0.18GeI/ 2 ,-B = —0.3GeIA, and a s is the strong coupling constant. The masses and 
parameters ( and £ for [c, q] and {c, q} are given in Table 1. S and A denote scalar and 
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axial vector diquarks of antisymmetric [c, q\ and symmetric {c, q} in flavor, respectively. 


The diquark is not indeed color-singlet, nor is it a real physical state. It generally 
exists in baryons. Thus, the estimated mass of a free diquark may differ from its actual 
mass in baryons. It is because of QCD interactions of the extra quark with the diquark 

pi 

Klciv et al. m discussed the uncertainty of mass prediction of charm-light diquarks 
as resulting from the uncertainties in QCD parameters. They presented the lower bound 
of [c, q] mass in two states 0 + and 1 + as 1.86±0.05 and 1.87±0.10, respectively, while they 
presented the upper bound of the mass as 2.02 and 2.07 for 0 + and 1 + states, respectively. 


Table 1: Masses and form factor parameters of heavy-light diquarks M- 


Quark content 

Diquark type 

M(MeV) 

£ (GeV) 

C (GeV 2 ) 

M 

S 

1973 

2.55 

0.63 

(c, q} 

A 

2036 

2.51 

1.05 

[c, s] 

S 

2091 

2.15 

1.05 

{c,s} 

A 

2158 

2.12 

0.99 


A number of theoretical approaches have been proposed to study heavy diquark masses. 
Examples include Bethe-Salpeter equation [18] , constituent diquark model [19], and rel¬ 
ativistic quark model based on a quasipotential approach in QCD [20]. It is worth men¬ 
tioning that diquark states in the nuclear matter are of masses larger than those states 
in the vacuum. Ref. [21] presents these mass uncertainties for light-flavor diquark states. 


Inserting the data given in Table 1 and diagonalizing the kernel and finding A = 1 
eigenvalue, we can calculate the binding energy. Using these results and the relation of 
mass and binding energy in Eq.(16), we obtain tetraquark system masses. 


M 


m 1 +m 2 + 


E h 


(16) 


where mi ,m 2 are diquark and antidiquark masses, respectively. In Table 2 and Table 3, 
binding energies and masses for charm tetraquarks are presented. The masses obtained 
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via this method turn out to be in such good agreement with the experimental and the¬ 
oretical masses that the observed errors are insignificant. In Fig.2 and Fig.3 we have 
compared our results with Maiani et al. and Ebert et al. 


Table 2: cqcq Masses and binding energies obtained via our method compared with theoretical 
predictions and possible experimental candidates in different J PC states 


Tetraquark 

j p C 

The calculated 

The calculated 

Mass in 

Mass in 

Exp 



E b (MeV) 

mass (MeV) 

ref [31 (25 ] 

ref [15] 


SS 

0++ 

-242.86 

3703.14 

3723 

3812 


(. SA+AS) 

V2 

1++ 

-147.32 

3861.68 

3872 

3871 

v r 3871.4 122; 
3875.2 [223 

(SA+AS) 

x/2 

1+- 

-255.2 

3744.84 

3754 

3871 


AA 

0++ 

-232.92 

3839.08 

3832 

3852 


AA 

1+- 

-185.8 

3886.2 

3882 

3890 


AA 

2++ 

-135.25 

3946.75 

3952 

3968 

v r 3943(231 
^-3914.31241 


Table 2: cscs Masses and binding energies predicted via our method compared with theoretical 

results in different J PC states 


Tetraquark 

J PC 

The calculated 

E b (MeV) 

The calculated 

mass (MeV) 

Mass in 

ref [15] 

SS 

0++ 

-128.1 

4053.9 

4051 

(SA + AS)/V 2 

1++ 

-127.29 

4111.71 

4113 

{SA-AS)/V 2 

1+- 

-150.34 

4098.66 

4113 

AA 

0++ 

-216.38 

4099.62 

4110 

AA 

1+- 

-180.47 

4135.53 

4143 

aA 

2++ 

-101.97 

4214.03 

4209 


It must be notified that we did not introduce any new or different parameters in potential 
for calculating the mass spectrum of heavy tetraquarks but employed the values already 
presented in [15| [18] and obtained the tetraquark mass spectrum by the two-body system’s 
binding energy. A good agreement was observed between our results and experimental 
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data and other references nasi 125]. Therefore, we can be confident that the potential 
coefficients chosen from ref. [15] in our calculations to solve Lippman-Schwinger equation 
were appropriate. Furthermore, we can conclude that our proposed binding energy was 
appropriate because it led to results nearly the same as those obtained by other scholars in 
previous studies. In Fig.2 and Fig.3, we have compared our results with other theoretical 
results for these systems. 


Figure 2: Mass spectrum of cqcq. a:our work, b:Ebert et al.’s work, c: Maiani et al.’s 
work 


In the recent work of Maiani et al. (type-II diquark model), diquarks are more 
resembling compact bosonic building blocks f2Bj. They have considered only diquark- 
antidiquark spin interactions and thus they have neglected spin-spin interactions between 
different diquarks. This means these results are an approximation of ref[3], Ignoring 
the spin interactions within diquarks and within antidiquarks and using different K- 
coefficients have contributed to the differences between the numerical value of the masses 
from experimental results and from their previous results [3]. 

The potential used in the paper is spin-independent. The mass differences observed 
in different J PC states result from the different masses and form factor parameters of 
heavy-light diquarks that we used for these states (Table 1). Thus, spin indirectly affects 
the tetraquark systems’ binding energies and masses (Table 2). Therefore, the mass 
differences observed between our results and those of ref. [15] derive from the fact that we 
ignored spin interaction in the potential. In our future works, we will add spin to the 
potential. 
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Figure 3: Mass spectrum of cscs. a: our work, b: Ebert et al.’s work 


5 Conclusion 

In this paper, we made use of the potential coefficients proposed by Ebert et al. and solved 
Lippman-Schwinger equation for heavy tetraquark systems. We managed to obtain the 
binding energy and used it to calculate the masses of heavy charm tetraquarks. The 
tetraquark is considered as the bound state of a heavy-light diquark and antidiquark. We 
used the diquark-antidiquark picture to reduce a complicated four-body problem to two 
simpler two-body problems. Our masses are very close to experimental and theoretical 
masses. Thus our method is appropriate for investigating heavy tetraquarks. Our results 
are in good agreement with the results derived from complicated relativistic methods and 
can be a good replacement for them. 

This method could equally be used for bottom tetraquarks and the tetraquarks composed 
of two heavy quarks and two heavy antiquarks. 
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